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METHOD AND APPARATUS FOR COINCIDENCE 
IMAGING DIGITAL TRIGGERING 

BACKGROUND OF THE INVENTION 

[0001] This invention relates generally to nuclear medicine imaging 
systems and, more particularly, to methods and apparatus for providing coincidence 
imaging in such systems. 

[0002] Positrons are positively charged electrons that are emitted by 
radionuclides that have been prepared using a cyclotron or other device. These 
positions are employed as radioactive tracers called "radiopharmaceuticals" by 
incorporating them into substances, such as glucose or carbon dioxide. The 
radiopharmaceuticals are injected into a patient and become involved in such 
processes as blood flow, glucose metabolism, and protein synthesis. 

[0003] Positrons are emitted as the radionuclides decay. The 
positrons normally travel a very short distance before they encounter an electron, and 
when this occurs, the positron and electron are annihilated and converted into a pair of 
photons, or gamma rays. This annihilation "event" is typically characterized by two 
features that are pertinent to positron emission tomography (PET) scanners, namely, 
each gamma ray has an energy of 5 1 1 keV and the pair of gamma rays are directed in 
substantially opposite directions. An image is created by determining the number of 
such annihilation events at each location within the field of view. 

[0004] A PET scanner may include two or more solid-state or 
scintillation detectors to detect individual photons. Some known scanners include a 
plurality of detectors that define a ring around a volume of interest. Timing the 
detection of these events is used to identify the pairs of photons from a single 
annihilation. To facilitate effective operation, detection events should be able to be 
timed to ten nanoseconds or less. Each scintillation detector includes a scintillator 
that converts the energy of each 511 keV photon into a flash of light that is sensed by 
a photomultiplier tube (PMT). Coincidence detection circuits are coupled to the 
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detectors and record only those photons that are detected simultaneously by two 
detectors located on opposite sides of some part of the patient. The number of such 
simultaneous events indicates the number of positron annihilations that occurred 
along a line joining the two opposing detectors. Within a few minutes hundreds of 
millions of events may be recorded to indicate the number of annihilations along lines 
joining pairs of detectors in the ring. These numbers of events are employed to 
reconstruct an image using, for example, computed tomography techniques. 

[0005] At least some known coincidence imaging systems use analog 
electronic circuits to generate a timing signal resulting from the sharp rise in output 
voltage in the detection device when the detection device detects an event. However, 
such circuits may have a limited accuracy when the event rate in the detector is very 
high. Typical trigger circuits often fail to detect events occurring immediately after a 
preceding event, and when they do detect such events, they detect an event time later 
than the actual time. 

BRIEF SUMMARY OF THE INVENTION 

[0006] In one aspect, a method for operating a coincidence imaging 
system is provided. The method includes receiving samples from a detector output 
and determining an intersection of a first line corresponding to a baseline portion of 
the detector output samples and a second line corresponding to a pulse rise portion of 
the detector output samples. 

[0007] In another aspect, a trigger circuit for a coincidence imaging 
system is provided. The circuit includes an analog-to-digital converter coupled to a 
detector output, a set of N accumulators coupled in a circuit parallel to the output of 
the analog-to-digital converter wherein N is a predetermined value, and a shift register 
coupled to the output of said plurality of accumulators. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0008] Figure 1 is a block diagram of an exemplary PET scanner 

system. 
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[0009] Figure 2 is a graph of an exemplary event detection pulse 
output of the PMT shown in Figure 1 . 

[0010] Figure 3 is an expanded view of the exemplary event 
detection pulse output shown in Figure 2. 

[0011] Figure 4 is a schematic diagram of an exemplary trigger 
circuit that may be used with the imaging system shown in Figure 1 . 

[0012] Figure 5 is a flowchart of an exemplary method that may be 
used to implement the trigger circuit shown in Figure 4. 

DETAILED DESCRIPTION OF THE INVENTION 

[0013] Figure 1 is a block diagram of an exemplary PET scanner 
system 5 including a gantry 1 0 which supports a detector ring assembly 1 1 about a 
central opening or bore 12. A patient to be examined is positioned in front of the 
gantry 10 and is aligned with the central axis of the bore 12. A motorized patient 
table (not shown) moves the patient into the bore 12 in response to commands 
received from an operator workstation 15. A gantry controller 17 responds to 
commands received from operator workstation 15 through a communication link 18, 
such as, for example, a serial communication link to operate gantry 10. 

[0014] In this exemplary embodiment, detector ring assembly 11 
comprises a plurality of radiation detectors 20. An exploded view of one radiation 
detector 20 illustrates that each radiation detector 20 includes a set of scintillator 
crystals 21 (for example, BGO crystals) arranged in a matrix and disposed in front of 
one or more photomultiplier tubes 22 (PMT) or other light detector. Each PMT 22 
produces an analog signal on at least one of a plurality of conductors 23 when a 
scintillation event occurs. A set of acquisition circuits 25 are mounted within gantry 
10 to receive these analog signals and produce digital signals indicating the event 
coordinates (x,y) and the total energy the scintillation event. The digital signals from 
each acquisition circuit 25 are sent through a cable 26 to a respective event locator 
circuit 27 that may be housed in a separate cabinet from gantry 10. Each acquisition 
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circuit 25 also produces a timing signal, which indicates when the scintillation event 
took place. 

[0015] The event locator circuits 27 form a portion of a data 
acquisition processor 30 that periodically processes the signals produced by 
acquisition circuits 25. Data acquisition processor 30 may have an acquisition CPU 
29 that controls communications on communication link 18 and a backplane bus 31. 
Event locator circuits 27 process and assemble the information regarding each valid 
event into a set of digital numbers (e.g., event data packet) that include a time marker 
that indicates when the event took place and the position of scintillation crystal 21 that 
detected the radiation event. Each event data packet is transmitted to a coincidence 
detector 32 that also forms a portion of data acquisition processor 30. 

[0016] Coincidence detector 32 receives the event data packets from 
event locator circuits 27 and determines if any two of the event data packets are in 
coincidence. Coincidence is determined by a number of factors. In one embodiment, 
the time markers in each event data packet must be within 12.5 nanoseconds of each 
other (or some other value of time window), and second, the locations indicated by 
the two event data packets must lie on a substantially straight line which passes 
through the field of view in scanner bore 12 and as described in more detail herein. 
Events that cannot be paired are discarded, but coincident event pairs satisfying the 
coincidence requirements are located and recorded as a coincidence data packet that is 
conveyed through a link 33 (e.g., serial link) to a sorter 34. In another embodiment, 
the timing signals of two events which are deemed to be in coincidence are used to 
calculate a time difference. The time difference is the time between the detection of 
the first event of a coincidence pair and the detection of the second event of the pair. 

[0017] Figure 2 is a graph of an exemplary event detection pulse 200 
output from PMT 22 (shown in Figure 1). Pulse 200 includes a baseline portion 202, 
a pulse rise portion 204, and a decay portion 206. An x-axis 208 represents time and 
a y-axis 209 represents an amplitude event detection pulse 200. An analog-to-digital 
converter (not shown) of acquisition circuits 25 digitizes the voltage output from PMT 
22 at a rate sufficient to obtain a plurality of samples 210, 212 and 214 during pulse 
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rise portion 204. In an exemplary embodiment, an event start 216 is detected by 
measuring a voltage difference between successive samples. Because noise may be a 
significant limitation in gamma ray detection, an event detection algorithm, executing 
on system 5, detects three successive voltage differences of greater than five percent 
of the full-scale voltage. In an alternative embodiment, other algorithms may 
selectably be used. 

[0018] In this exemplary embodiment, a temporal accuracy of less 
than one nanosecond may be achieved by determining the time at which a straight line 
218 fitted through pulse rise portion 204 intersects a straight line 220 fitted through 
baseline portion 202. 

[0019] Figure 3 is an expanded view of the exemplary event 
detection pulse 200 (shown in Figure 2). A difficult event detection situation occurs 
during relatively high event rates when a new event starts before the previous event is 
completed. In an exemplary embodiment, an accurate estimate of the actual event 
time (e.g., the time the gamma ray entered detector 20) is obtained by determining the 
time where a straight line 302 fitted through baseline portion 202, prior to start of the 
event, intersects a straight line 304 fitted through pulse rise portion 204. 

[0020] Pulse rise portion 204 may be fitted to the equation, 

y = m x t + b u 

where m x is the slope and 

b\ is the y-intercept. 

[0021] In the exemplary embodiment, a y-axis 306 (t = 0) (shown as 
a vertical line) arbitrarily may be chosen on pulse rise portion 204, and slope m } at 
that point may be determined as the average difference between the measurements 
immediately preceding and following the chosen y-axis. In an alternative 
embodiment, a time at which an event detector (not shown), using successive 
differences, detects the start of an event is selected as y-axis 306. Intercept, b u with 
this arbitrary y-axis is the average value of points immediately prior to and 
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immediately following, y-axis 306. In this exemplary embodiment, the baseline slope 
may be measured similarly, for example using y = m 2 t + b 2 , where m 2 is the slope or 
the baseline, and b 2 is the y-intercept.. In an alternative embodiment, the baseline 
slope may be assumed to be zero. The baseline intercept, b 2 , is then the average value 
of measurements made prior to the start of the rise portion 204. 

[0022] The time of the start of the pulse is then calculated from: 
X = -(b\-b 2 )/m x 

[0023] Result, t, is the time between the start of the event (the point 
where lines 302 and 304 intersect) and the arbitrary y-axis 306, in units of the 
sampling interval. This method may be embodied within a software code segment 
executing on acquisition CPU 29 within system 5 or may be implemented using 
dedicated firmware or hardware. 

[0024] Figure 4 is a schematic diagram of an exemplary trigger 
circuit 400 that can be used with imaging system 5 (shown in Figure 1). Trigger 
circuit 400 includes an analog-to-digital converter 402 including an input 404 and an 
output 406. Input 404 is coupled to an output of a photon detector, such as PMT 22. 
Output 406 is coupled to an input of a plurality of accumulators 408 that are coupled 
in parallel. An output of each accumulator 408 is coupled to an input of a shift 
register 410 that includes a plurality of register elements 412. Shift register 410 
includes a plurality of outputs, Ql, Q2, Q3, and Q4 that are selected to generate 
proper delay times for computational accuracy of the pulse samples taken during 
baseline portion 202 and rise portion 204 of pulse 200. Ql represents the sum of N 
values immediately prior to arbitrary y-axis 306. Q2 represents the sum of N values 
immediately after arbitrary y-axis 306. Q3 represents the sum of N values centered 
on arbitrary y-axis 306. Q4 represents the sum of N values prior to the start of rise 
portion 204. 

[0025] Ql and Q2 are transmitted to inputs of a computational logic 
circuit, such as an adder circuit 414 that combines Ql and Q2 into an output, and in 
particular, the subtraction of Ql from Q2 (Q2-Q1). Q3 and Q4 are transmitted to a 
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computational logic circuit, such as an adder circuit 416 that combines Q3 and Q4 
into an output Q3-Q4. An output of adder circuit 414 is coupled to an input of a look- 
up table component 418. An output of look-up table component 418 and an output of 
the adder 416 are combined by a multiplier 420. 

[0026] In operation, analog-to-digital converter 402 transmits voltage 
measurements at a substantially constant rate to the set of N accumulators 408 and 
shift register 410. Accumulators 408 each sum N successive voltage measurements 
and output the sums to shift register 410, such that successive elements 412 of shift 
register 410 contain the following summed voltage measurements: 

N N N N N 

'-' 1=1 '=1 i=l i=l i=| 

For example, in one embodiment where N = 4, 
suml = vl+v2+v3+v4 
sum2 = v2+v3+v4+v5 
sum3 = v3+v4+v5+v6 
sum4 = v4+v5+v6+v7 

[0027] When an event (i.e., coincidence event) is detected, for 
example, by a threshold method, four terms, Ql, Q2, Q3, and Q4 are extracted from 
shift register 410 for processing. A delay in time from coincidence measurement to 
extraction at the Q4 element is determined by a response of a threshold detector (not 
shown) and is selected to be sufficient to ensure that in all circumstances, Q4 contains 
only values prior to rise portion 204. Q4, therefore is the sum of N samples from the 
ADC prior to the start of the pulse rise, so represents the baseline immediately prior to 
the start of the rise. Q3 is the sum of N samples centered on some arbitrary time, T, 
during the rise, so the average of these values should equal the signal value at this 
arbitrary time. Q2 is the sum of N samples some time after the arbitrary time, T, and 
Ql is the sum of N samples preceding the arbitrary time T by the same duration as Ql 
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lags T. If N were chosen to be 1, then Q3 - Q4 would be the signal rise from the 
baseline to time T, and Q2 - Ql would be the rise from the time Ql was collected to 
the time that Q2 was collected. Dividing Q2-Q1 by the time from Ql to Q2 gives the 
slope of the rising part of the pulse. 

[0028] The value of Q3-Q4 represents the value N*(b,-b 2 ) and the 
value of Q2-Q1 represents N values of (V i+k -Vj) where k is the number of sampling 
intervals between the start of the Ql and the Q2 sums. If N is selected as a power of 2 
(e.g., 2, 4, 6, 8, etc.) then a factor of N can be removed by shifting or by selecting 
only the most significant bits for subsequent processing. The resulting values Bl = 
(b, - b 2 ) = (Q3-Q4)/N and Ml = (Q2-Q1)/N are the intercept and slope of the straight 
line fitting the samples Ql to Q3, where the y values are relative to the average 
baseline, b 2 . The intersection of the rising portion of the pulse with the baseline is 
the value of (Q3-Q4)/(Q2-Q1) is determined using by look-up table 418 and 
multiplier 420. This is the time, relative to the arbitrarily defined, T, at which the 
rising slope intercepts the baseline. 

[0029] In an exemplary embodiment, using a sodium iodide 
scintillation detector yields an event rise time of approximately eighty to ninety 
nanoseconds and a decay time of about eight hundred nanoseconds (with a two 
hundred thirty nanosecond half-time). Analog-to-digital converter 402, sampling at 
200 MHz (i.e., every five nanoseconds) generates approximately sixteen sample 
values during event pulse rise portion 204. The event detector (not shown) examines 
differences in every fourth sample, and triggers if three successive such differences 
exceed a predetermined limit, such as five percent of full scale. Such a limit provides 
a trigger time about thirty to forty nanoseconds after the start of the event, which 
becomes the arbitrary y-axis 306 (shown in Figure 3). 

[0030] The time calculation uses N=8. With an 8-bit analog-to- 
digital converter 402, each of the Qn outputs has a range often bits, assuming, for 
example, that the first half of the sampled rise is within the lower half of full scale. 
The difference Q2-Q1 yields nine bits. The most significant eight bits are used to 
access a look-up table component 418 with a ten bit inverse value (i.e., one bit whole 
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and nine fractional bits), which in turn is multiplied, using multiplier 420, by the ten 
bit value, Q3-Q4. The resulting twenty-bit value contains eleven integral bits and 
nine fractional bits. However, this result is N times the true value oft, and so, if N=8, 
the result should be interpreted as fourteen integral and six fractional bits. This time 
value may be entered directly as the event time. 

[0031] A precise time marker may be generated at a time exactly T d 
after the start of the pulse. In the exemplary embodiment, the arbitrary y-axis 306 
(shown in Figure 3) is about 30-40 nanoseconds after the start of the event, 
accordingly T d may be selected to equal, eighty nanoseconds (or sixteen sample 
intervals). Time value, t, calculated as described above, is the time from the start of 
the event to the arbitrary y-axis 306, thus, the timing pulse should be generated at 16-t 
intervals after the arbitrary y-axis. (16-t) is calculated and the whole number part 
loaded into a countdown counter that counts down using the sample clock and 
generates a timing pulse at zero. This pulse (timed to the preceding clock before the 
desired time) is then provided to a programmable delay line, which has been loaded 
with the fractional part of (16-t), such that the pulse is delayed by the fractional part. 
The integral part of the result can be used to control a countdown timer that is 
partitioned into three whole number bits, and three or more fractional bits. 

[0032] Figure 5 is a block diagram of an exemplary method 500 that 
may be used to implement trigger circuit 400 (shown in Figure 4). Method 500 
includes receiving 502 samples from a detector output and determining 504 an 
intersection of a first line corresponding to a baseline portion 202 of the detector 
output samples and a second line corresponding to a pulse rise portion 204 of the 
detector output samples. 

[0033] Exemplary embodiments of apparatus and methods that 
facilitate timing coincidence events in a PET scanner are described above in detail. A 
technical effect of the timing apparatus and methods described herein include at least 
one of facilitating improving determination of coincidence pairs of positron 
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annihilation events. The timing apparatus and methods include a plurality of logic 
devices and a processor that allows for controlling a timing trigger and logic outputs. 

[0034] It will be recognized that although the system in the disclosed 
embodiments comprises programmed hardware, for example, executed in software by 
a computer or processor-based control system, it may take other forms, including 
hardwired hardware configurations, hardware manufactured in integrated circuit form, 
firmware, among others. It should be understood that the trigger circuit disclosed 
may be embodied in a digital system with periodically sampled signals, or may be 
embodied in an analog system with continuous signals, or a combination of digital 
and analog systems. 

[0035] The above-described methods and apparatus provide a cost- 
effective and reliable means for facilitating generating a time signal to perform 
coincidence imaging. More specifically, the methods and apparatus facilitate 
improving the accuracy of time signals during periods of relatively high count-rates. 
As a result, the methods and apparatus described herein facilitate operating 
coincidence imaging systems in a cost-effective and reliable manner. 

[0036] Exemplary embodiments of coincidence imaging systems are 
described above in detail. However, the systems are not limited to the specific 
embodiments described herein, but rather, components of each system may be utilized 
independently and separately from other components described herein. Each system 
component can also be used in combination with other system components. 

[0037] While the invention has been described in terms of various 
specific embodiments, those skilled in the art will recognize that the invention can be 
practiced with modification within the spirit and scope of the claims. 
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